Group analysis and exact solutions of a class 

of variable coefficient nonlinear telegraph equations 

Ding-jiang Huang and Nataliya M. Ivanova ^ 

^ Department of Applied Mathematics, Dalian University of Technology, Dalian, 
116024, P- R- China 
e-mail: hdj8116@163.com 

Institute of Mathematics of NAS of Ukraine, 3 Tereshchenkivska Str., 01601 Kyiv, Ukraine 
e-mail: ivanova@imath.kiev.ua 

A complete group classification of a class of variable coefficient (l+l)-dinicnsional telegraph 
equations f{x)utt — {H{u)ux)x + Kiu)ux, is given, by using a compatibility method and 
additional equivalence transformations. A number of new interesting nonlinear invariant 
models which have non-trivial invariance algebras are obtained. Furthermore, the possible 
additional equivalence transformations between equations from the class under consideration 
are investigated. Exact solutions of special forms of these equations are also constructed via 
classical Lie method and generalized conditional transformations. Local conservation laws 
with characteristics of order of the class under consideration are classified with respect to 
the group of equivalence transformations. 

1 Introduction 

Since Sophus Lie (1842-1899) introduced the notion of continuous transformation group, now 
known as Lie group, the theory of Lie groups and Lie algebras have been evolved into one 
of the most explosive development of mathematics and physics throughout the past century. 
Nowadays, this theory has been widely apphed to diverse fields of mathematics including dif- 
ferential geometry, algebraic topology, bifurcation theory, numerical analysis, special functions 
and to nearly any area of theoretical physics, in particular classical, continuum and quantum 
mechanics, fluid dynamics system, relativity, and particle physics [7,10,18,20-22,35,37,50]. 

When applied to system of differential equations, the mathematical trends whose object is a 
common treatment of the Lie groups of transformations and the differential equations admitted 
by these groups is called group theoretical analysis of differential equation. Traditionally, there 
are two interrelated problems in this subject. The first one is finding the maximal Lie (sym- 
metry) transformation group admitted by a given equation. The second problem is classifying 
differential equations that admit a prescribed symmetry group G. The principal tool for han- 
dling both problems is the classical infinitesimal routine developed by S. Lie (see, e.g., [35,37]). 
It reduces the problem to finding the corresponding Lie symmetry algebra of infinitesimal oper- 
ators whose coefficients are found as solutions of some over-determined system of linear partial 
differential equations. 

The problems of group classification and exhaustive solutions of such problems are not only 
interesting from the purely mathematical point of view, but also important for applications. 
It is well known that modelling the phenomena in nature (such as in physics, chemistry and 
biology) by partial differential equations is one of the central problems of mathematical physics 
and applied mathematics. Generally, those modelling differential equations could contain some 
arbitrary parameters or functions which have been found experimentally and so are not strictly 
fixed. Therefore, in order to reflect the natural laws accurately one has to decide which dif- 
ferential equation fits in the best way as a model for the process under study and so has to 
select from a broad class of possible partial differential equations. Because in many physical 



models there often exist a priori requirements for symmetry groups that follow from physical 
laws (in particular, from Galilean or relativistic theory), which imply that solving the problems 
of group classification makes it possible to accept for the criterion of applicability the following 
statement: modelling differential equations have to admit a group with certain properties or 
the most extensive symmetry group from the possible ones. This point of view is supported by 
the fact that the most successful mathematical models in theoretical and applied science have a 
rich symmetry structure. Indeed, the basic equations of modern physics, the wave, Schrodinger, 
Dirac and Maxwell equations are distinguished from the whole set of partial differential equa- 
tions by their Lie and non-Lie (hidden) symmetries (see, e.g., [18] for more details on symmetry 
properties of these equations). 

In the approach used here, an exhaustive consideration of the problem of group classification 
for a parametric class L of systems of differential equations includes the following steps: 

1. Finding the group Cx^'^'^ (the kernel of maximal Lie invariance groups) of local transforma- 
tions that are symmetries for all systems from L. 

2. Construction of the group C (the equivalence group) of local transformations which 
transform L into itself. 

3. Description of all possible G~-inequivalent values of parameters that admit maximal in- 
variance groups wider than Cp-^^ . 

Following S. Lie, one usually considers infinitesimal transformations instead of finite ones. This 
approach essentially simplifies the problem of group classification, reducing it to problems for 
Lie algebras of vector fields. See [1,2,34,37,41,55] for precise formulation of group classification 
problems and more details on the used methods. 

The result of application of the above algorithm is a list of equations with their Lie invariance 
algebras. The problem of group classification is assumed to be completely solved if 

i) the list contains all the possible inequivalent cases of extensions; 

ii) all the equations from the list are mutually inequivalent with respect to the transformations 
from C^; 

Hi) the obtained algebras are the maximal invariance algebras of the respective equations. 

Such a list may include equations that are mutually equivalent with respect to local trans- 
formations which do not belong to Knowing such additional equivalences allows one to es- 
sentially simplify further investigation of £. Constructing them can be considered as the fourth 
step of the algorithm of group classification. Then, the above enumeration of requirements for 
the resulting list of classifications can be completed by the following step: 

iv) all the possible additional equivalences between the listed equations are constructed in 
explicit form. 

In this paper we consider a class of variable coefficient (l-l-l)-dimensional nonlinear telegraph 
equations of the form 

j{x)uu = {H{u)Ux)x + K{u)ux (1) 

where / = f{x), H = H{u) and K = K{u) are arbitrary and sufficient smooth real- valued 
function of their corresponding variable, f{x)H(u) ^ 0. In what follows, we assume that 
(Hu,Ku) 7^ (0,0), i.e., (1) is a nonlinear equation. This is because the linear case of (1) 
{H, K = const) has been studied by Lie [30] in his classification of linear second-order PDEs 
with two variables. (See also a modern treatment of this subject in [37]). 
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The study of equation (1) is stimulated not only their intrinsic theoretical interest, but also its 
physical importance. Equations (1) are used to model a wide variety of phenomena in physics, 
chemistry, mathematical biology etc. For the case f{x) = 1 and K(u) = equation (1) can be 
used to describe the flow of one-dimensional gas, longitudinal wave propagation on a moving 
threadline and dynamics of a finite nonlinear string and so on [3,4]. When K{u) = this 
equation describes the longitudinal vibrations of an elastic and non-homogeneous taut string or 
bar [52]. The outstanding representative of the class of equations (1) is the nonlinear telegraph 
equation that is the mathematical model for a large number of physical phenomena. (For more 
details refer to [3,31].) 

Historically, there are a number of papers contributed to the studies of Lie groups of trans- 
formations of various class of (1 -|- l)-dimensional nonlinear wave equations and their individual 
members. Probably, Barone et al [6] was the first study of the following nonlinear wave equation 

utt = Uxx + F{u), 

by means of symmetry method, this equation was also studied by Kumei [28] and Pucci et al [45] 
subsequently. Motivated by a number of physical problems, Ames et al [3,4] investigated group 
properties of quasi-linear hyperbolic equations of the form 

Utt = [f{u)ux]x- (2) 

Later, their investigation was generalized in [16, 23, 52] to equations of the following forms 
respectively 

Utt = [f{x, u)ux]x, Utt = [fiu)ua: + g{x, u)]a:, and uu = f{x, u.^)uxx + 9{x, u^). 

The alternative form of equation (2) was also investigated by Oron and Rosenau [36] and Suhubi 
and Bakkaloglu [51]. Arrigo [5] classified the equations 

Uu = U^U:cx + F{u). 

Furthermore, classification results for the equation 

Utt + K{u)ut = [F{u)U:c\x 

can be found in [21]. An expand form of the latter equation 

Utt + K{u)ut = [F{u)Ua:]x + H{u)Ux 

was studied by Kingston and Sophocleous [27]. Recently, Lahno, Zhdanov and Magda [29] 
presented the most extensive list of symmetries of the equations 

Utt = Uxx + F{t, X, u, Ux) 

by using the infinitesimal Lie method, the technique of equivalence transformations and the 
theory of classification of abstract low-dimensional Lie algebras. There are also some papers [15, 
19, 36, 44] devoted to the group classification of the equation of the following form 

Utt = F{uxx), Utt = F{ux)uxx + H{ux), and utt + Xu^x = 9{u, u^). 
It worthwhile mentioned that the equations 
Utt = {F{u)ux)x + H{u)ux 
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together with its equivalent potential systems have also been studied by Bluman et al [8,11-13]. 
In their a series of papers, many interesting results (especially for case of power nonlinearities) 
including Lie point and nonlocal symmetries classification and conservation laws of the four 
equivalent systems were systematically investigated. 

Prom the above introduction, it is easy to see that equation (1) is different from any aforemen- 
tioned ones. However, equation (1) is a generalization of many physically important systems, 
thus there is essential interest in investigating them from a unified and group theoretical point. 
The ultimate goal of this paper is to present an extended group analysis and to find additional 
equivalence transformations and exact solutions of equations (1). A lot of new interesting cases 
of extensions of the maximal Lie symmetry group and cases with high- dimensional spaces of 
conservation laws were obtained for these equations. 

Problems of general group classification, except for really trivial cases, are very difficult. 
This can be illustrated by the multitude of papers where such a general classification problem 
is solved incorrectly or incompletely. There are also many papers on "preliminary group classi- 
fication" where authors list some cases with new symmetry but do not claim that the general 
classification problem is solved completely. Por this reason, finding an effective approach to 
simplification is essentially equivalent to showing the feasibility of solving the problem at all. 
Recently, based on the investigation of the specific compatibility of classifying conditions, Nikitin 
and Popovych [34] developed an effective tool (we refer it as compatibility method) for solving 
the group classification problem of nonlinear Schrodinger equation. Their method has been 
applied to investigating a number of different group classification problem [14,25,34,40,41,53]. 
In particular, in [41] Popovych and Ivanova extended the method to complete group classifica- 
tion of nonlinear diffusion-convection equations by further considering the so called additional 
equivalence transformations. However, to the best our knowledge, there are no any result about 
application the compatibility method to equation (1). Therefore, the paper is one of new appli- 
cation of the compatibility method to the problem of group classification. Although the authors' 
debt to the works of Nikitin, Popovych and Ivanova [34,41] is evident, the results of group clas- 
sification of class (1) presented in this work seem to be new. Hence, these will lead to some 
explicit applications in Physics and Engineer. 

The rest of this paper is organized as follows. Since the case f{x) = 1 has a great variety 
of applications and has been investigated earlier by a number of authors, we collect results for 
this class together in Section 2. In Section 3 we present the results of the complete group clas- 
sification of class (1). Some additional equivalence transformations are considered in Section 4, 
where we also present the result of group classification of class (1) with respect to the set of 
point transformations. The result of the group classification is used to find exact solutions of 
equations from class (1) (Section 5). Ibid we construct functionally separation solutions for 
some equations (1). A natural continuation of group analysis of equations (1) is investigation of 
their conservation laws. More precisely, in Section 6 we construct the local conservation laws of 
equations of form (1) having characteristics of order 0. Finally, some conclusion and discussion 
are given in Section 7. 

2 Group classification for the subclass with f{x) = 1 

Class (1) includes a subclass of equations of the general form 

utt = {H{u)Ux)x + K{u)ux (3) 

(i.e., the function / is assumed to be equal to 1 identically). These equations (called the 
nonlinear telegraph equations) are important for applications. Symmetry properties of class (3) 
were studied by a number of authors (see [13,19,21,27] for details). However, in all the above 
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references the results of group classification of class (3) are presented in implicit form only. 
Therefore we single out the results of the group classification of equations (3) from classification 
of class (1). 

Theorem 1. The Lie algebra of the kernel of principal groups of (3) is A\*^^ = {dt, dx)- 
Theorem 2. The Lie algebra of the equivalence group for class (3) is 

Ai = {dt, dx, du, udu, tdt - 2HdH - 2KdK, xdx + 2HdH + Kdx). 
Any transformation from GJ" for the class (3) is 

t = te4 + ei, x = xe5 + e2, u = ueQ + e^, H = He'l^e\, K = Ke^'^e^, 
where ei, . . . , £7 are arbitrary constants, £46^6^ 7^ 0. 

Theorem 3. The complete set of -inequivalent extensions of 7^ A'^*^'" for equation (3) 
is exhausted by ones given in table 1. 

Table 1. Case of f(x) = 1 



N 


H{u) 


K{u) 


Basis of A™^ 


1 


V 


V 


dt, dx 


2 


V 





dt, dx, tdt + xdx 


3 




e"" 


dt, dx, (p. - 2i')tdt + 2(/x — v)xdx + 'idu 


4 


e" 


1 


dt,dx, tdt + 2xdx + 2du 


5 







dt, dx, tdt - 2du, xdx + 2du 


6 


\< 




dt, dx, {pL — 2v)tdt + — v)xdx + 2udu 


7 




1 


dt, dx, l^tdt + 2nxdx + 2udu 


8a 


\< 





dt, dx, jJbtdt — 2udu, nxdx + 2ud„ 


8b 






dt, dx, tdt + udu, e~''{dx + udu) 


9 


-4 

U 


-4 

u 


dt, dx, 2tdt + udu, t^dt + tudu 


10 







dt, dx, 2tdt + udu, 2xdx — udu, t^dt + tudu 


11 







dt, dx, 2tdt + 3udu, 2xdx — 3udu, x^dx — 3xudu 



Here {fi, v) ^ (0, 0), {fj,, v) € {(1, 0), (0, 1)} mod G~ in case 3, fj, ^ 0, —4, —4/3 in case 8. Case 8b is reduced to 
case 8a by means of transformation i = t, x = , u = e~^u. 

Remark 1. The most similar form of classification results have been presented in [13]. Note 
that the point symmetries of some special classes of equation (3) adduced in table III of [13] can 
be obtained, respectively, by appropriately scaling t and u or taking some appropriately values 
for the parameters n and u in cases 1.3 and 1.6. For example, the classes A and E of table III 
in [13] can be obtained by setting /x = q;+1,i^ = 1 and n = a, u = 1 respectively in case 1.3. 

Remark 2. The proof of theorem 3 follows directly from the results of the next section. 

3 Results of Classification 

Following the above-mentioned algorithm we are looking for an infinitesimal operator in the 
form 

Q = T{t, X, u)dt + ^{t, X, u)dx + r]{t, x, u)du (4) 

which corresponds to a one-parameter Lie group of local transformations and keeps the equa- 
tion (1) invariant. The classical infinitesimal Lie invariance criterion for equation (1) to be 
invariant with respect to the operator (4) read as 

pr(2) Q(A) |a=o= 0, A = f{x)uu - {H{u)ux)x - K{u)ux. (5) 
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Here pr(^) Q is the usual second order prolongation [35,37] of the operator (4). Substituting the 
coefficients of pr(^) Q into (5) yields the following determining equations for r, ^ and rj: 



^t = iu = Vuu = 0, Hr]^^ + Krjx - fvu = 0, 



ifH^ - 2TtH - r,Hu + 2Hi^ = 0, 

i^Cxx - 2H^j]^ - j]Ku - 2TtK + ^K^ + ^^K - 27]^uH = 0, 
2vtuf - Tttf = 0, 2HC^u - 2TtHu + ^Hu^ + 2F„^^ - r/„F„ - r/i7„„ = 0. 

Investigating the compatibility of system (6) we find that the final equation of system (6) is an 
identity (substituting the second equation of system (6) to the final one can yield this conclusion). 
With this condition, system (6) can be rewritten in the form 

Tx = Tu=it=^u= Vuu = 0, 2ritu = Tu, (7) 

2(e. - n) + = (8) 

Hrjxx + Kijx - fr]tt = 0, (9) 

{H^K - K^H) j^-K^x- 2H^Vx + H{^^^ - 2r/^„) = 0. (10) 

Equations (7) do not contain arbitrary elements. Integration of them yields 

r = T{t), ^ = C(x), r] = r]\t,x)u + r]\t,x), 7j\t , x) = + aix) . (11) 

Thus, group classification of (1) reduces to solving classifying conditions (8)-(10). 

Splitting system (8)-(10) with respect to the arbitrary elements and their non-vanishing 
derivatives gives the equations = 0, = 0, 77 = for the coefficients of the operators from A^^^ 
of (1). As a result, the following theorem is true. 

Theorem 4. The Lie algebra of the kernel of principal groups of (1) is A^^'^ = (df). 

The next step of the algorithm of group classification is finding equivalence transformations 
of class (1). An equivalence transformation is a nondegenerate change of the variables t, x and u 
taking any equation of the form (1) into an equation of the same form, generally speaking, with 
different /(x), H{u) and K{u). The set of all equivalence transformations forms the equivalence 
group C . To find the connected component of the unity of C, we have to investigate Lie 
symmetries of the system that consists of equation (1) and some additional conditions, that is 
to say we must seek for an operator of the Lie algebra A~ of G~ in the form 

X = Tdt + + vdu + -Kdf + pdH + A (12) 

from the invariance criterion of (1) applied to the system: 

f{x)utt = {H{u)ux)x + K{u)ux, 

ft = fu = 0, Ht = H^ = 0, Kt = K^ = 0. ^^^^ 

Here u, f, H and K are considered as differential variables: u on the space {t, x) and /, H, G on 
the extended space (t, ,t, u). The coordinates r, r/ of the operator (12) are sought as functions 
of i, X, u while the coordinates tt, p, 99 arc sought as functions of t, x, u, f, H, K. 

The invariance criterion of system (13) yields the following determining equations for r, ^, ry, 
TT, p and ip: 

Tx = Tu = = ^u = Vx = Vu = 0, Ttt = Vuu = 0, 

TT* = TTu = TTfl- = TTK, Pt = Px = Pu = Pf = PK, ^t = Vx = = 0, 

j + 2i^- 2Tt -PH = 0, (14) 
^xxH+{j+^x-2n)K-ip = 0. 
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After easy calculations we find from (14) 

r = ci + c^t, ^ = C2 + C5X, ■>] = C3 + cqu, 

TT = (C7 + 2C4 - 2C5)/, p = C7H, (p= {C7- C5)K, 

where ci, . . . , C7 are arbitrary constants. Thus, we obtain the following statement. 

Theorem 5. The Lie algebra of for class (1) is 

A- = {dt, d^,du, tdt + 2fdf, xd^ - 2fdf - kOk, ud^, fOf + Hd^ + kOk). 

Continuous equivalence transformations of class (1) are generated by the operators from A". 

For class (1) there also exists a non-trivial group of discrete equivalence transformations gen- 
erated by four involutive transformations of alternating sign in the sets {t},{x,K}, {u} and 
{f,H,K}. Therefore, C contains the following continuous transformations: 

i=te4^ + ei, x = X€5 + e2, u = u€Q + e3, f = fe\e^^e-j, H = He-j, K = Ke-je^^, 

where ei, . . . , ey are arbitrary constants. 

Theorem 6. A complete set of inequivalent equations (1) with respect to the transformations 
from G"^ with A^^ 7^ A^^'^ is exhausted by cases given in tables 2-4- 



Table 2. Case of yH{u) 



N 


K(u) 


/(^) 


Basis of A-"*^ 


1 


V 


V 


dt 


2a 


V 




dt, etdt + 2d^ 


2b 


H 


g-2x-7e-"^ 


dt, 7t9t + 2e^a^ 


2c 


H 


e (e +7) 


dt, (A + 2)fat -2(l+7e")a, 


2d 





\x\' 


dt, (A + 2)tdt + 2x8^ 


3a 





1 


dt, d^, tdt +xdx 


3b 


H 




dt, tdt — dx, e^dx 



Here 7, A 0, e € {0, 1} mod G~, 7 = ±1 mod G~. 
Additional equivalence transformations: 

1. 2b 2&{K = 0, e = 1): t^t,x = -76""', u = u; 

2. 2c(A ^ -2) ^ 2a(i4: = -H/{\ + 2), e = 1): f = S = (A + 2) In I7 + e"''!, w = w; 
2c(A = -2) 2&{K = -H, e = 0): i = t, x = \ii\-f + e-"], u = u; 

3. 2d(A 7^ -2) ^ 2a.{K = -H/{X + 2), e = 1): t = t,x = {X + 2) In |a;|, ti = u; 
2d(A = -2) ^ 2a(is: = -/f, e = 0): i = t, x = \n\x\, u = u; 

4. 3b ^ 3a: i = t, x = e~^, u = u. 



Table 3. Case of H{u) = e" 



N 


A* 


K(u) 




Basis of A-""" 


1 


V 






dt, 


[A(m -u) + {fi- 2u)]tdt + 2{n - v)xdx + 2du 


2 


V 




1 


dt. 


dx, {p — 2v)tdt + 2(/i - v)xdx + 2du 


3 


1 


1 




dt, 


(A + l)tdt + 2x9:, + 2du 


4 


1 


1 


1 


dt. 


dx, tdt + 2xdx + 2du 


5 


1 


ee" 


V 


dt. 


tdt - 2du 


6a 


1 





/'(^) 


dt, 


tdt - 2du, atdt + 2{l3x'^ + 71a; + 70)^;^ + 2l3xdu 


6b 


1 


e" 




dt. 


tdt - 2du, atdt - 2(7oe" + 71 + /3e-")9^ + 2/3e-"a„ 


7a 


1 





1 


dt, 


tdt — 2du, xdx + 2du, dx 


7b 


1 


e" 




dt. 


tdt - 2du, dx - 2du, e'dx 


7c 


1 







dt. 


tdt — 2du, xdx — du, x^dx + xdu 


7d 


1 


e" 


e-2-(e— +7)-3 


dt. 


tdt - 2du, tdt + 2(1 + 7e")a., (e-" + -ife^'dx - (e"" + 7)9„ 
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Here (ju, i^) # (0,0), (/^, z/) € {(1, 0), (0, 1)} mod G~; 7 = ±1 mod G~; A # 0; e € {0,1} mod G~; a, /3, 71, 
70 = const and 

/3e-^ - a - 27oe"= 



/ (a;) = exp 



-3/3x - 271 + a 



da; 



/ (a;) = exp 



Px'^ + 'iix + 70 
Additional equivalence transformations: 

1. 6b 6a: t = t, x 

2. 7b ^ 7a: i = t, x 

3. 7c ^ 7a t = t signx, i = 1/a;, = w — In 

4. 7d^7a: f=tsign(7 + e-^), x = 1/(7 + 6"'^), M = M-ln|(7 + e-^)| 

Table 4. Case of H{u) = lul" 



/Je-'^ +71+70 



-dx 



N 








Basis of A™"" 


1 


V 


\u\'' 


\x\^ 


a, [A(^ - i^) + (/i - 2!/)]fat + 2(^ - i^)xd^ + 2wa„ 


2 


V 


\< 


1 


9t, dx, - 2u)tdt + - u)xdx + 2udu 


3 


V 


1 




dt, (A + \)iitdt + 2/ia;(?a; + 2vBu 


4 


V 


1 


1 


dt, dx, ntdt + 2iixda: + 2udu 


5 


#-4 




V 


dt, (itdt — 2udu 


6a 


#-4 







dt, ntdt — 2udu, atdt + 2[{fi + l)Px^ + 71a; + 7o]5x + 2/3xudu 


6b 


#-4 






dt, )j,tdt - 2udu, 










atdt - 2[(Ai + l)^e-" + 71 + 7oe"]S^ + 2j3e-''udu 


7a 


#-4,-1 





1 


dt, ntdt — 2udu, dx, nxdx + 2udu 


7b 


-4 -4 






dt, iitdt — 2udu, tidx — 2udu, e'dx 


7c 


# -4,-1,-1 

' ' 3 ' 





|a;| f+i 


dt, fJ,tdt — 2udu, ti{lJ. + V}xdx — (m + 2)vBu, 










+ l)x^dx + xudu 


7d 


_4 _4 _i 




— 2a: 

e 


dt, ntdt - 2udu, ifi + 2)tdt + 2{n + l)(e-" + -t)e'dx, 




















(At + l)(e-" + -yfe^'dx - (e"" + -i)udu 


7e 


-1 







dt, tdt + 2udu, dx — judu, tdt + xdx — jxudu 


7f 


-1 






dt, tdt + 2udu, e'dx — judu, tdt — dx + ■ye~''udu 


7g 


-2 




1 


dt, dx, tdt + udu, e~''{dx + udu) 


8 


-4 




V 


dt, 2tdt + udu, t^dt + tudu 


9a 


-4 





/'WU=-4 


dt, 2tdt + udu, t^dt + tudu. 










atdt + 2[— 3/3a;^ + 71a: + 7o]9a, + 2Pxudu 


9b 


-4 






dt, 2tdt + udu, t^dt + tudu. 










atdt - 2[-3/3e-" + 71 + 7oe"]a, + 2/3e-"Ma„ 


10a 


-4 





1 


dt, 2tdt + udu, t^dt + tudu, dx, 2xdx - udu 


10b 


-4 






dt, 2tdt + udu, t^dt + tudu, 2dx + udu, e'^dx 


10c 


-4 





8 

a; 3 


dt, 2tdt + udu, t^dt + tudu, 6xdx + vBu, Zx'^dx — xudu 


lOd 


-4 






dt, 2tdt + udu, t^dt + tudu, tdt + 3(e-^ + f)e''dx, 








(e-»=+7)3 










3(e-+7)2e-a. + (e-^+7)ua„ 


11a 


4 
3 





1 


9t, 2tft + 3udu, dx, 2xdx — 3udu, x^dx — Sxudu 


lib 


4 
3 


4 

M 3 




dt, 2tdt + 3udu, 2dx + 3udu, e~^{dx + 3udu), e^dx 



Here (fx, i^) # (0, 0), A / 0, e € {0, 1} mod G"; # for cases 3-7d; a, Co, /3, 71, 70 = const and 

-(3/i + 4)/3x - 271 + a^^\ ^4,„, _ / f {H + 2)Pe-'' -a- 27oe" 

(/u + l)/3e-^ + 71 + 7oe^ 



f{x)=expi^j ^j^'^ ' /*(»:) = exp -i / T '. 7sT „ . . dx 



l)f3x^ + 71a; + 70 
Additional equivalence transformations: 

1. 6b 6a, 7b^7a, 9b 9a, 10b 10a, lib ^ 11a : i=t,x = e""=, u = u: 

2. 7c 7a, 10c 10a (// = -4): i=t,x = -i, -u = \x\~^u. 

3. 7d ^ 7a, lOd 10a (a* = -4): i=t,x = it = I7 + e-^'pTTT^w. 

4. 7e — > 7a (/i = —1): t = f, 5 = x, -u = e^'^u; 

5. 7f ^ 7a (m = -1): i = t,x = e"' , u = e-^'=""w. 

6. 7g ^ 7a (/i = -2): f = t, x = e^, w = e-'^M. 
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Proof. To prove the theorem we use the compatibility method [34,41]. The basic idea of this 
method is based on the fact that the substitution of the coefficients of any operator from 
^max y ^kcr jj^^^ ^j^g classifying equations results in nonidentity equations for arbitrary elements 
(see [34,41] for more details and exhaustive examples of applications). In our case the procedure 
of looking for the possible cases mostly depends on equation (8). For any symmetry operator 
equation (8) gives some equations on H of the general form (an + b)Hu = cH, where a, b, c 
are constant. For all operators from A™'^ the number k of such independent equations is not 
greater than 2; otherwise they form an incompatible system on H. k is an invariant value for 
the transformations from C. Therefore, there exist three inequivalent cases for the value of k: 
i) A; = 0: H{u) is arbitrary; ii) k = 1: H{u) = or H{u) = ^ 0) mod and iii) k = 2: 

H(u) = 1 mod Let us consider in more detail case H{u) = (table 3). We attempted to 
present our calculations in reasonable detail so that verification would be feasible. For this case 
equations (8) and (11) imply r)u = 0, i.e. r) = i]{t,x) and tu = 0. Therefore, equations (8)-(10) 
can be written as 

2(?x - Tt) + J^ = V, (15) 

e"r?,, + Kr,^ - frfu = 0, (16) 
{K - Ku)v - KCx - ^e^Vx + e^ixx = 0. (17) 

Equation (17) looks like Ku = + 6e" with respect to K, where u,b = const. Therefore, K 
must take one of the following four values. 

(i) K = e*^" + /iie" mod where i^G{0, 1}, hi = const. Substituting K into equations (16) 
and (17) yields rj = const, hi = and = {l^ — i^)v which implies there exist two cases for /, 
i.e. / 7^ const or not by further considering equation (15). Thus we get cases 1 and 2. 

(ii) K = ne" + /lie" mod G~ and hi = const. It follows from equations (16) and (17) that 
T] = for any operator from A^^^, which contradict with the assumption that f] ^ 0. 

(iii) K = + ho mod C, where ho = const. Substituting K into (16) and (17) yields 

Vxx + Vx = 0, horjx - fr]tt = 0, Cxx - Cx - 2% = 0, ho{r] - Cx) = 0. (18) 

Solving the first and the third equation of system (18) we obtain r] = C^(t)e~^ + C'^(^)) ^ = 
7oe^ + 7i — C^{t)e^^. Since = 0, we have = 0. Then it follows from the second and the 
fourth equation of system (18) and equation (15) that ho = and 

T = ^{C2 - a)t + ci, ^ = 7oe^ + 7i+/3e-^, r/ = -/^e"^ - C2, 

where ci, C2, a, /?, 70, 71 = const. Hence, equation (15) implies that the function / must satisfy 
Z (/ = 0, 1, 2) equations of the form 

fx Pe-"" -a- 27oe^ 



/ /3e =^ + 71 + 7oe^ 



with non-proportional sets of constant parameters (a, /3, 70, 71). The values I = and I = 1 
correspond to cases 5 (e = 1) and 6b. An additional extension of yl™^^ exists for I = 2 m 
comparison with Z = 1 iff / is a solution of the equation 

fx A2e-^ 



/ Aie-^ + Ao ' 

where either A2 = or A2 = 3Ai 7^ 0. Integrating the latter equation gives cases 7b and 7d. 
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(iv) K = ho = const mod In an analogous way to that in the previous case, we obtain 
ri = C^{t)x + C°(i)) C = li^ + 70 + C^{t)x'^, what is more, ^, rj satisfy 

hoVx - fvtt = 0, ho{r] - = 0. 
Investigating the compatibiUty of the latter system and equation (15) with = leads to 

T = ^{c2 + a)t + ci, ^ = + 7ia; + 7o, rj = Px - C2, 

where ci, C2, a, /?, 70, 71 = const, and (3, C2, 71 satisfy ho^ = 0, ho(c2 + 71) = 0. Hence, there 
exist two cases for ho, i.e. /iq = or not. The vahic /iq = 1 mod G~ results in cases 3 and 4. 
Below. h(] = 0. Equation (15) holds when the function / is a solution of a system of l{l = 0, 1, 2) 

equations of the form 

_ -3/?x + g - 271 

/ + 71X + 70 

with non-proportional sets of constant parameters (a, /3, 70, 71). The values I = and I = 1 
correspond to cases 5 (e = 0) and 6a. Additional extension of exists for Z = 2 in comparison 
with Z = 1 if and only if / is a solution of the equation 

fx _ ^2 
/ Aix + Ao 

where either A2 = or A2 = — 3Ai / 0. These possibilities result in cases 7a and 7c. 

The rest of the cases of values H can be studied in an analogous way. □ 

In what follows, for convenience we use double numeration T.N of classification cases and 
local equivalence transformations, where T denotes the number of the table and N the number 
of the case (or transformation) in table T. The notation 'equation T.N' is used for the equation 
of the form (1) where the parameter functions take the values from the corresponding case. 

The operators from tables 2-4 form bases of the maximal invariancc algebras if the cor- 
responding sets of the functions /, H, K are G~-inequivalent to ones with most extensive 
invariance algebras. For example, in case 4.1 (//,z^) / (0,0) and A 7^ —6 \i u = 1. Similarly, in 
case 3.1 the constraint set on the parameters /x, u and A coincides with the one for case 4.1, and 
|x = 1 if 1/ = 0. 

4 Additional equivalence transformations and classification 
with respect to the set of point transformations 

In tables 2-4 we list all possible C-inequivalent sets of functions /(x), H{u), K{u) and corre- 
sponding invariance algebras. However, these tables contain some cases being equivalent with 
respect to point transformations that do not belong to The simplest way to find such addi- 
tional equivalences between previously classified equations is based on the fact that equivalent 
equations have equivalent maximal Lie invariance algebras. 

Explicit formulas for additional transformations that do not change the value of H{u) are 
adduced after the tables. Besides these transformations there exist additional point transforma- 
tions changing H{u). Thus, e.g., 

t = X, X = t, u = Inu 

maps case 4.7a to 3.7a. One more example of similar transformations is 

i=x, x = t, u = u'^'^^, fi = —fi/{fi + 1) 
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between equations of form utt = {u^^Ux)x-, A* 7^ ~1- I^i particular, it connects cases 4.10a 
and 4.11a. The same transformation applied is a discrete symmetry for equation with ^ = —2. 
The latter two transformations are, indeed, partial cases of more general transformation 

f=.. i = t, * = (19) 

between equations from class 

uu = {H{u)ux)x, (20) 

where the new transformed value of arbitrary element H is the derivative to the inverse function 
u = H{u) for = / H{u)du [21]. Transformation (19) is nonlocal with respect to the arbitrary 
element H{u) and therefore can be considered as generalized extended equivalence transforma- 
tion [24,32] in class of nonlinear wave equations (20). 

One can check that there exist no other point transformations between the equations from 
tables 2-4. Using this we can formulate the following theorem. 

Theorem 7. Up to point transformations, a complete list of extensions of the maximal Lie 
invariance algebra of equations from class (1) is exhausted by the cases given in table 2, cases 
3. 1-3. 6a and 4- 1-4- 10a numbered with Arabic numbers without Roman letters and subcases "a" 
of each multi-case. (Two equations from case 4-7 with parameter values /x and fx are assumed 
to be equivalent iff ft = — /// (^ + 1) )■ 

As one can sec, the above additional equivalence transformations have multifarious structure. 
This displays a complexity of a structure of the set of admissible transformations. Usually the 
problems of finding of all possible admissible transformations are very difficult to solve, see, 
e.g., [26,27,39,43]. We will try to discuss the structure of the set of admissible transformations 
of class (1) in a sequel paper. 

A more systematic way to proceed with equivalence transformations is to classify them using 
the infinitesimal method or the direct method. Examples of conditional equivalence algebras 
calculated by the infinitesimal method are listed in table 5. 

Table 5. Conditional equivalence algebras 



Conditions 



Basis of A" 



K = H 
K = H = e^ 
H = e", K = 
K = H = u'' 
H = u'',K = 



dt,d^, tdt + 2/9/, du, udu, e"(9, - 2fdf), fdf + Hdn 

dt, d^, tdt + 2/9/, du + fdf, c"(a. - 2fdj), e-^-d^ +du- fdf) 

dt, dx, tdt + 2/9/, du + fdf, xdx - 2/9/, x'^dx + xdu — ^xfdf 

dt, dx, tdt + 2fdf, udu+iifdf, e'{dx-2fdf), e-"[(l + At)9x -m9„ + (2 + //)/9/] 

dt, dx, tdt + 2/9/, udu + nfdf, xdx - 2/9/, (1 + n)x'^dx + xudu - (4 + i^i)xfdf 



To find the complete collection of additional local equivalence transformations including both 
continuous and discrete ones, we should use the direct method. Moreover, application of this 
method allows us to describe all the local transformations that are possible for pairs of equations 
from the class under consideration. A problem of this sort was first investigated for wave 
equations by Kingston and Sophocleous [26,27,49]. 

Remark 3. It is a well-known that class (20) is linearizable [9] with respect to potential hodo- 
graph transformation (e.g., interchange of independent and dependent variables) applied to the 
potential system Vx = ut, vt = H(u)ux corresponding to the simplest local conservation law 
of (20) . This transformation can be considered as potential equivalence transformation between 
the classes of wave equations uu = F{x)uxx and (20). 
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5 Exact solutions 



In this section, we turn to the presentation of some exact solutions for (1) by means of the 
classical Lie-Ovsiannikov algorithm and generalized conditional symmetry methods. We first 
present the solutions of some special forms of the nonlinear wave equation (20). Then using 
our classification with respect to all the possible local transformations, we transform them to 
solutions of more complicated telegraph equations (such as 3.7b, 4.7b, 4.7e). 



5.1 Exact solutions obtained via clcissical Lie-Ovsiannikov algorithm 

Let us note that the equations with / = 1 are well investigated and that most of the exact 
solutions given below have been constructed before (see citations in [13,19,21,27]). However, to 
the best of our knowledge, there exist no works containing a systematic study of all the possible 
Lie reductions in this class, as well as exhaustive consideration of the integrability and exact 
solutions of the corresponding reduced equations. That is why we have decided to implement 
the relevant Lie reduction algorithm independently, especially since it is not a difficult problem. 
So, let us consider equation 3.7a: 

utt = {e^Ux)x- (21) 

Let us recall that for (21) the basis of A"^^^ is formed by the operators 

Qi = dt, Q2 = tdt - du, Q3 = dx, = xdx + 2du. 

The only non-zero commutators of these operators are [Qi, (52] = Qi and [Qs, Q4] = Qs- There- 
fore is a realization of the algebra 2^2.1 [33]. All the possible inequivalent (with respect 
to inner automorphisms) one-dimensional subalgebras of 2^2. 1 [38] are exhausted by the ones 
listed in table 6 along with the corresponding ansatzc and the reduced ODEs. 

Table 6. Reduced ODEs for equation (21). a ^ 0, e = ±1. 



N 


Subalgebra 


Ansatz u = 


U! 


Reduced ODE 


1 


(Qi> 




X 


(e"^)" = 


2 


(Q2) 


(fi{io) - 21n[i| 


X 


(e*')" = 2 


3 


(Q3) 




t 


<p" = 


4 


(Q4> 


ip{ui) + 2ln\x\ 


t 


(e"^)" = 2e'^ 


5 


(Qi + eQa) 




X — et 


(e^)" = eV" 


6 


(Q2+eQ3) 


ip{io) - 2ln\t\ 


X — eln \t\ 


(e^)" = eV" + evj' + 2 


7 


(Qi + eQi) 


if>{uj) + 2et 




(e^)" = + v?') 


8 


{Q2 + aQi) 


ip{Lj)+2{a - l)ln|t| 


x\t\- 


(e'^)" = a^^^V" + a(a + l)W - 2(a - 1) 



We succeeded in solving the equations 6.1-6.5. Thus we have the following solutions of (21): 



u = In Icix -|- co|, u = In 



X + ClX + Co 



U = Cit + Co, 



u 



= '"(4cgcosh\^^)+"')' - = ^(--^^)' 



where if satisfies e'^ = e'^ip + ci(p + cq. Using these we can construct solutions for cases 3.7b- 
3.7d easily. For example, the transformation 3.4 yields the corresponding solutions for the more 
complicated and interesting equation (case 3 .7d) 



e ^""{e + 7) ^utt = {e!'u^)x + e^u^ 



(22) 
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in the following form 

■u = In I cix + co(e~^ + 7)1) u = ln 



+ S + S(e-^ + 7) 



U = Cit + Cq. 



The power /x = — 1 is a singular value of the parameter /j, for case 4.7a. So, the corresponding 
equation 



Utt = {u ^Ux) 



(23) 



is distinguished by the reduction procedure. It is remarkable that cases 4.7e and 4.7f are reduced 
exactly to equation (23). Exact solutions of equation (23) can be easily obtained by direct 
application of the classical Lie reduction method or using transformation (19) applied to solutions 
of equation (21). Both these approaches lead us to the following solutions of (23): 



u = C2e'^i^, u = — 



u = x {t"^ + C2t + ci) , u = (p{x — et), 



4cf cosh2(i^; 



U = Clt + Co, u 



4c2cos2(^)' 



where satisfies J ]r-^dip = co + ci. Analogously to the previous case, we obtain by means of 



C2'P 

transformations 4.5 exact solutions of equation 4.7f, i.e. 
in the following forms: 



(24) 



u = C2e 



(ci +7)e 



U 



4c2cOSh2(^^ 



U = Clt + Co, u 



4c2cos2(^^tC2- 



u 



-2x+'ye 



where is as above. 

Another example of a variable coefficient equation is given by case 4.7d. To look for exact 
solutions of it, first we reduce it to equation 

Utt = {u^Ux)x (25) 
(case 4.7a). As in the previous cases, the invariance algebra of (25) is of the form 

Qi = dt, Q2=tdt-2^r^udu, Q3 = dx, Qa = xdx + 2ir^udu. 

It is also a realization of the algebra 2A2.i- The reduced ODEs for equation (25) are listed in 
table 7. 

Table 7. Reduced ODEs for equation (25). a / 0, e = ±1. 



N 


Subalgebra 


Ansatz u — 


UJ 


Reduced ODE 


1 


(Qi) 


LP{U}) 


X 




2 


(Q2> 




X 


W)' = f(f + 1V 


3 


(Qs) 


(p{uj) 


t 


(^" = 


4 


(Q4) 


(p{oj)\x\'i^ 


t 




5 


(Qi+eQs) 


ip{ui) 


X — et 




6 


(Q2 + eQs) 


(p{uj)\t\'~ 


X — e In \t\ 




7 


(Qi + eQ4,} 






(vjV')' = e't^V" + 4£'m"V - (4/X-' - l)e^u>y:>' 


8 


{Q2 + aQi) 


2(c-l) 

(p{uj)\t\ 


x\t\- 


{ip'^^'Y = a^t^V" - a[(4a - 4)^-^ - a- l]wip' 
+2{a - l)/i-i[(2Q - 2)/x-i - l]v5 
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For some of the reduced equations we can construct the general solutions. For others we 
succeeded in finding only particular solutions. These solutions are the following: 

1 

u = |c2 — cix — ci/ix| 1+'' , u = cit + co, u = ip{x — et), 

where ip satisfies j^^^^'^ — (p — ciui — C2 = 0. All the results of table 7 as well as the solutions 
constructed can be extended to equations 4.7b-4.7f using the local equivalence transformations. 
So for the equation (case 4.7d) 

e ^^(e ^ + 7) M+i utt = {u^u^)x + u^'u^ 
the transformations 4.3 yield exact solutions in the form 

-u = |c2(e~^ + 7) + ci + ci/x|iT^, u = cit + co, u = \{e~'' + -f)\T^ip{-{e'^ + jy'^ - et), 
with the same value of tp. 

5.2 Functionally separation solutions obtained via generalized conditional 

symmetry methods 

Wc now turn to the functionally separation solutions of the nonlinear telegraph equations (3) by 
using the generalized conditional symmetry methods [17,46,54]. This method was developed by 
Fokas, Zhdanov and Qu et al, and has been applied to study the functional separation of variables 
for various nonlinear equation [47, 48] . In order to implement the method effectively, let us 
review some basic notations of the functionally separation solutions and generalized conditional 
symmetry [46-48]. 

Definition 1. A solution u{t, x) of equation (3) is said to be functionally separable if there exist 
functions q{u), ip{t), and ^(x) such that 

qiu)=^{t)+ijix), (26) 

where q{u) is some smooth function of u, (pit) and ip{x) are some undetermined functions of t 

and X respectively. 

The classical additively separable solution and product separable solution are particular cases 
of the above functional separable solution. 

Definition 2. An evolutionary vector field 

V = r]{t,x,u,...)du 

is said to be a generalized conditional symmetry of (3) if 

wV^^\uu - {H{u)u^)^ - K{u)u^)\Enw = 0, (27) 

where E is the solution manifold of (3) and is a second-order system of (3) obtained by ap- 
pending the condition 77 = and its differential consequences; pr F^^) is the second prolongation 
of the infinitesimal operator V. 

According to reference [47], there exist the following theorem. 

Theorem 8. Equation (3) possesses the functional separable solution (26) iff it admits the 
generalized conditional symmetry 

q"iu) 

V = {u:ct + 9{u)uxUt)du, g = -7^- (28) 

q [u) 
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Substituting (28) into (27) and using (3) , a straightforward calculation shows that the func- 
tions H, K and g satisfy the following system 

H"-gH'^^ = Q, K'-K^^;- = 0, g" - 2gg' + {g^ - g')M^ = 0, 

{2gg' - g")H + {g^ - g')H' - 2gH" - 2g'H' + H'" + gH" 

-{gH' + H")^+H'i3g'-2g^) = 0, 

K" - g'K - {gK + K')^ + K{y - 2g^) = 0. 

To obtain solutions of this system, we consider three cases for the third equation. 

(A) g' — g^ = 0. In this case, g is given hj g = and g = —1/u after by translating u. 
Substituting g = into the remain equations of system (29) and scaling u, we find that 

H = e"", K = ce"", 

where a, c is arbitrary constant. We always assume a 7^ to exclude the cases equation (3) is 
linear. Furthermore, from (28), we have q = u. 
For g = —1/u, we can obtain 

H = u°', K = cu'^, q = lnu. 

(B) g' — g^ / 0, g" — 2gg' = 0. In this case, it is easy to see that H, K = const. Equation (3) 
is linear. We thus do not consider it. 

(C) g' — g^ 7^ 0, g" — 2gg' ^ 0. In order to solve system (29), we define h{u) such that 
g = —h'/h. From the first and the third equation of system (29), we obtain 

where a is nonzero constant (if a = equation (3) becomes linear) , Hq is a constant and can be 
chosen as ±1 by scaling t; h satisfies hh'" — h'h" = ah", which can be integrated as 

hh" - = ah' + b. (31) 

Since a 7^ 0, there exist two cases for solving equation (31), i.e., 6 = or not. If 6 = 0, we can 
find from (30), (31) and the rest equations of system (29) that 

' ^-"a/(dc) + e'^«' 

where d is arbitrary constant. By scaling and translating u, we can set d = ±1 and a/c = ±1. 
Hence, four possibilities are distinguished: 

/i = c(e-"-l), H={e''-l)-\ K = 0, g = ln(e" - 1), 

when Hq = —d = 1, a/c = 1; 
/i = c(e" + l), = e"(e" + l)-\ K = 0, g = -u - ln(e" + 1), 

when Hq = —d = —1, a/c = 1; 
/i = c(e"-l), i7 = e"(e" - l)-\ K = 0, g = -n - ln(e" - 1), 

when Hq = —d = —1, a/c = —1; 
/i = c(e-« + l), i/=(e" + l)-\ K = 0, q = ln{e'' + l), 

when Hq = —d = —1, a/c = —1. 
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For 6 7^ we introduce F{h) = in (31), F satisfies an ordinary differential equation 
hFFfi = F'^ + aF + b, which is integrated to 



zd:. , / h 



In 



^— ^, /lo = const. (32) 



z'^ + az + b \ho 

Set A = — 46, we must consider three cases separately for the above integral. 
For A = 0, we have 

which can be rewritten as F = h' = pi{h). 
For A > 0, we have 

+ VA^ Va+« _|_ o - VA^ Va-« - ^ ^ j 

We rewrite it implicity as h' = p2{h). 
For A < 0, we have 

, /2/i' + a\2] r ^ a /2/i' + o\i //i\2 

1 + — , exp — 2 , arctan — , = i— ■ 

We rewrite it implicity as h' = ps (h) . 

In these cases, h, H can be determined implicity by 

dz pj + api + b . 
^ = u, H{u) = , z = 1,2,3, 

and K satisfies the second and the fifth equations of system (29). The corresponding equation 

has a separable solution of the form 



^ dz 



Summing up the above analysis, we have the following results: 

Theorem 9. Equation (3) possesses the functional separable solution (26) if and only if it is 
G~ equivalent to one of the following equations: 

(a) utt = {e^Ux)x + ee"u^, u = (p{t) + tp{x), 

(b) Utt = {u^Ux)x + eu'^Ux, u = (p(t)il}{x), 

(c) = ((e" - 1)- V)., ln(e"-l) = (/.(i)+V(x), 

(d) Utt = (e"(e" + ir^ux)x, u - ln(e" + 1) = ip{t) + V(x), 

(e) utt = {e^{e^-l)-'ux)x, -u - \n{e^ - 1)) = ip{t) + ^{x), 

(f) n« = ((e" + 1)- V)., ln(e" + l)=^(i)+V(x), 

(g) Utt = i- ^' ^ ^ Ux)x + K(u)ux, J = (yf (t) + V'(a;), i = l,2,3, 
where e = 0, 1, fi ^ 0, pi are described above. 
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Using this theorem and additional equivalence transformations, we can construct some exact 
solutions for equations (3) and (1). For example, equation (a) in theorem 9 admits the functional 
separation solution u = ip{t) + ip{x), where tp and ^ satisfy the system 

cp" = Xe'^, (e'^VO' + ee'^^' - A = 0, 
which can be solved explicitly by 

ci 



<^(t) = In l^sec ^^^t + C2^^^ ' = InjAx + C3e ^ - C4 - A|, e = 1. 

If we take e = 0, the equation (a) is exactly Eq. (21) which admits the separation solution 



u = hi 



2Ac? 



+ In 



A 2 

-X - C3X + C4 



This solution is novel which can not be obtained classical Lie-Ovsiannikov algorithm. By using 
the transformation 3.4, we obtain an interesting functional separation solution of the equa- 
tion (22) as follows: 



= In 



sec(23^isign(7 + e '') + ^ 



C2 \2 



2Acf 



+ ln 



A 



2(7 + e-^) 



— :r-C3 + C4(7 + e ^) 



Equation (b) in theorem 9 admits the separation solution u = ip{t)ip(x), where if and 
satisfy the system 

^" = A^'^+i , (V'^V'O' + eV'^'V'' - AV' = 0. 
When taking = — 1, e = 0, equation (b) become (23) which admits a new separation solution 
1 /I2 C2\2/A2 \ 

Using transformation 4.5, we can get a solution for Eq. (24) 



2Ac2 



sec 



V2ci 2ci 



C2_\2 /A 
^c^J V2 



r + est + C4 . 



Similarly, equation (c) admits the functionally separation solution u = ln(l + if{t)i/j{x)), ip and 
ip satisfy the system 

ip'^ = aif^ - 2X<p - P, V'^ - /^V'^ + 2AV'^ + = 0. 
The implicit solution of this system is given by 

Mt) dz . /•'^(^) dy 



I 



^Jaz^ -2Xz- P 



= X. 



Equations (d), (e), (f) also admit the functionally separation solution u = —\n{(p{t)ip(x) — 1), 
— ln(l — ip{t)^{x)), ln(ip{t)ip{x) — 1) respectively, (p and ip satisfy the system 



„'2 



ip-~ = a(p^ - 2X(p - P, V + Pip + '^M - OiW = 0. 
The implicit solution of this system is given by 

Mt) dz f^^'^'^ dy 



f 



^Jaz^ -2Xz-p 



^/-Py^ - 2A2/3 + ay2 



X. 
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6 Conservation laws 



In this section we classify local conservation laws of equations (1) with characteristics depending, 
at mostly, on t, x and u. For classification we use the direct method described in [42] . 

To begin with, we adduce a necessary theoretical background on conservation laws, following, 
e.g., [35,42] and considering for simplicity the case of two independent variables t and x. See 
the above references for the general case. 

Let £ be a system L{t, x, U(^p^) = of Z PDEs = 0, . . . , L' = for the unknown functions 
u = {u^, . . . , u^) of the independent variables t and x. Here U(^p-j denotes the set of all the partial 
derivatives of the functions u of order not greater than p, including u as the derivatives of the 
zero order. 

First we give an empiric definition of conservation laws. A conservation law of the system C 
is a divergence expression 

DtT{t, x, U(^r)) + DxX{t, X, = (33) 

which vanishes for all solutions of C. Here Dt and arc the operators of total differentiation 
with respect to t and x, respectively. The differential functions T and X are correspondingly 
called a density and a flux of the conservation law and the tuple (T, X) is a conserved vector of 
the conservation law. 

The crucial notion of the theory of conservation laws is one of equivalence and triviality of 
conservation laws. 

Definition 3. Two conserved vectors {T,X) and (T',X') are equivalent if there exist func- 
tions T, X and H olt, x and derivatives of u such that T and X vanish for all solutions of C and 
T' = T + T + DxH, X' = X + X — DtH. A conserved vector is called trivial if it is equivalent 
to the zeroth vector. 

The notion of linear dependence of conserved vectors is introduced in a similar way. Namely, 
a set of conserved vectors is linearly dependent iff a linear combination of them is a trivial 
conserved vector. 

Conservation laws can be investigated in the above empiric framework. However, for deeper 
understanding of the problem and absolutely correct calculations a more rigorous definition of 
conservation laws should be used. 

For any system C of differential equations the set CV{C) of conserved vectors of its con- 
servation laws is a linear space, and the subset CVo(i2) of trivial conserved vectors is a linear 
subspace in CV(>C). The factor space CL(>C) = CV(>C)/ CVo(>C) coincides with the set of equiv- 
alence classes of CV(>C) with respect to the equivalence relation adduced in definition 3. 

Definition 4. The elements of CL(>C) are called conservation laws of the system jC, and the 
whole factor space CL(>C) is called the space of conservation laws of C. 

That is why description of the set of conservation laws can be assumed as finding CL(£) 
that is equivalent to construction of either a basis if dimCL(£) < oo or a system of generatrices 
in the infinite dimensional case. The elements of CV(£) which belong to the same equivalence 
class giving a conservation law are considered all as conserved vectors of this conservation law, 
and we will additionally identify elements from CL{C) with their representatives in CV(>C). For 
(T, X) G CV{jC) and T £ CL(£) the notation (T, X) G JT will denote that (T, X) is a conserved 
vector corresponding to the conservation law J^. In contrast to the order ?'(t,X) of ^ conserved 
vector (T, X) as the maximal order of derivatives explicitly appearing in the differential functions 
T and X, the order of the conservation law T is called min{r(7^ x) I (^)^) ^ •^}- Under linear 
dependence of conservation laws we understand linear dependence of them as elements of CL(>C). 
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Therefore, in the framework of "representative" approach conservation laws of a system L are 
considered as linearly dependent if there exists Hnear combination of their representatives, which 
is a trivial conserved vector. 

Let the system L be totally nondegenerate [35] . Then application of the Hadamard lemma to 
the definition of conservation law and integrating by parts imply that the left hand side of any 
conservation law of £ can be always presented up to the equivalence relation as a linear combi- 
nation of left hand sides of independent equations from L with coefficients A'' being functions 
oit^x and derivatives of u: 

DtT + D^X = \^L^ + --- + \^LK (34) 

Definition 5. Formula (34) and the /-tuple A = (A^, . . . , A') are called the characteristic form 
and the characteristic of the conservation law DtT + D^X = correspondingly. 

The characteristic A is trivial if it vanishes for all solutions of C. Since C is nondegenerate, 
the characteristics A and A satisfy (34) for the same conserved vector (T, X) and, therefore, 
are called equivalent iff A — A is a trivial characteristic. Similarly to conserved vectors, the 
set Ch(>C) of characteristics corresponding to conservation laws of the system £ is a linear 
space, and the subset Cho(£) of trivial characteristics is a linear subspace in Ch(£). The factor 
space Chf(i2) = Ch(i2)/ Cho(i2) coincides with the set of equivalence classes of Ch(i2) with 
respect to the above characteristic equivalence relation. 

An important property of the class of equations in the conserved form is that it is preserved 
under any point transformation (see, e.g., [42]). 

Proposition 1. A point transformation g: i = t^{t,x,u), x = x^{t,x,u), u = u^{t,x,u) 
prolonged to derivatives of u transforms the equation DtT + D^X = to the equation 
DfT^ + DxX^ = 0. The transformed conserved vector (T^ ,X^) is determined by the formula 

^a/~ ~ ^ T{x,U(^r))Dti+ X{x,U^r))Dxi 
T^(x,U(r\) = = — , 

^ ^"^'^ DttD^x - DJDtx 

T{x, U(^j.))DtX + X{x, U(^r))DxX 



X^{x,U^r)) 



Dft DxX — Dxt DfX 



Remark 4. In the case of one dependent variable (m = 1) g can be a contact transformation: 
t = t^{t,x,U(i^), X = x^{t,x,U(i^)^ = it^^^(t, re, tt(i)). Similar note is true for the below 
statement. 

Proposition 2. Any point transformation g between systems C and C induces a linear one- 
to-one mapping g^ from CY{jC) into CV(£), which maps CVo(>C) into CVo(>C) and generates a 
linear one-to-one mapping gt from CL{C) into CL(>C). 

In such way, if a point transformation connects two systems of partial differential equations, 
then the same transformation maps the set of conservation laws of the first system to the set 
of conservation laws of the second system and the space of characteristics of conservation laws 
to the space of characteristics. Therefore, a group of equivalence transformations of a class of 
systems establishes a one-to-one correspondence between conservation laws of systems from the 
given class. So, we can consider a problem of investigation of conservation laws with respect 
to the equivalence group of a class of systems of differential equations This problem can be 
investigated in the way that is similar to group classification in classes of systems of differential 
equations. Namely, we construct firstly the conservation laws that are defined for all values 
of the arbitrary elements. (The corresponding conserved vectors may depend on the arbitrary 
elements.) Then we classify, with respect to the equivalence group, arbitrary elements for each 
of that the system admits additional conservation laws. 
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For more detail and rigorous proof of the correctness of the above definitions and statements 
see [42]. 

Using the most direct method described in [42] we prove the following theorem. 

Theorem 10. A complete list of -inequivalent equations (1) having nontrivial conservation 
laws with characteristics of the zeroth order is exhausted by ones given in table 8. 



Tabic 8. Conservation laws of equations (1) 



N 


H{u) 


K{u) 




Basis conservation laws 


1 


V 


V 


V 




2 


V 





V 


CL\ CL^, CL3, CL* 


3 


V 


H 


V 


CL\ CL2, CL^, 


4 


V 


1 


1 


CL\ CL^ CL^ CL* 


5 


V 


1 




CL\ CL^, CL^, CLi° 


6 


V 


H + k 




clS CL^, CL", CL12 


7 


V 


H-e 


(l + ce-")-i 


CLS CL^, CLl^ CL^^ 


8 


V 


H + P 


(l + ce-")-i 


CL\ CL^ CLl^ CL^'^ 



Here c and k are arbitrary constants, fc 7^ in cases 7 and 8. 

The conserved densities T and fluxes X of the presented conservation laws have the following 
forms: 

CL^: fut,-{Hu^ + jKy, 

CL^: f{tut-u),-t{Hu^ + jKy, 

CL^: xfut, -{xHua; - J H); 

CL^: x/(itxt - u), -t{xHu^ - J H); 

CL^: e^/ut, -e'^Hu^; 

CL^: e^/(t«t - u), -eHHu^; 

CV: {2x -t^)ut + 2tu, -{2x-t^){Hu^ + u) + 2 J H; 

CL^: (6te - - (6x - 3t2)tt, -(6te - t^){Hu:c +u)+6tjH; 

CL'^: sin tut — cos tu, — sint{xHux + xu — J H); 

CL^'^: cos tut + sin tu, — cos t{xHux + xu — J H); 

CLii; (2e2^ - kt^e'')ut + 2kte'^u, -(2e2^ - kt^e''){Hux + A;u) + H^e^ / H- 

CLi2. (6te2=^ - kt^e^)ut - i&e^x - Zkfe^)u, -(Gte^^ - kt^e^){Hux + ku) + kt^e'' J H; 

CLi3. e^-^*(«t + few), -e-^*(e^ + c){Hux +u)~ ce'''* J H; 

QLi4. ^x+kt^^^ _ j^^-j^ _e^'*(e^ + c){Hux + u) - ce^^ J H; 

CL^^: e^ {sin ktut — k cos ktu) , — sin kt{e^ + c)(Hux + u) — c sin kt j H; 

CL^^: (cos fetiit + sin /sill) , — cos kt{e^ + c){Hux + u) — c cos kt J H. 

The above conservation laws can be used for construction of potential systems, potential 
symmetries and potential conservation laws. We will present such analysis elsewhere. 

7 Conclusion 

In summary, we have performed completely group classification of the class of equations (1) by 
using the compatibility method and additional equivalence transformations [34,41]. The main 
results on classification are collected in tables 2-4 where we list inequivalent cases of extensions 
with the corresponding Lie invariance algebras. Following the tables, we write down all the ad- 
ditional equivalence transformations, reducing some equations from our classification to others 
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of simpler forms. For a number of equations from the list of the reduced ones we construct 
optimal systems of inequivalent subalgebras, corresponding Lie ansatze and exact solutions. By 
means of additional equivalence transformations the solutions obtained are transformed to the 
ones for the more interesting and complicated variable coefficient equations. Functionally sepa- 
ration solutions are obtained for a number of equations from class (1) via generalized conditional 
symmetry method. 

The present paper should be an inspiration for further investigations of different properties of 
class (1). For example, one can classify the nonclassical (conditional) symmetries. Furthermore, 
one can solve the general equivalence problem for any pair of equations from class (1) with 
respect to the local transformations so as to finding the group of all possible local equivalence 
transformations (i.e. not only continuous ones) in the whole class (1) as well as all the conditional 
equivalence transformations. 

Motivated by some elegant results about potential symmetries of partial differential equations 
by Bluman et.al [8-13] and results of section 6, we also intend to systematically calculate nonlocal 
symmetries and higher order local and potential conservation laws, construct invariant and 
nonclassical solutions, as well as obtain linearizations, etc, by investigating nonlocally related 
potential systems and subsystems of variable coefficient equations (1). Further work along these 
lines would be extremely interesting. 
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